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We consider a mean-field Hamiltonian for a 0^,2_y 2 -|- (p -|- ip) superconductor(SC) in presence 
of spin-density-wave(SDW) order. This is due to the fact that the non-commutativity of any two 
orders produces the third one. The energy spectrum of such a Hamiltonian is shown to be gapped 
and it yields a topological phase in addition to the conventional one. A phase diagram characterizing 
different topological phases is construted. The Chern numbers and hence the nature of the topo¬ 
logical phases are determined. The edge state spectrum and the possibility of whether the vortex 
state harbouring the zero modes are discussed. 


The subject of topological systems in condensed mat¬ 
ter is one of the most active field of research at present 
and is developing in rapid pace [HQ. Topological phases 
are characterized by the existence of both the gapless 
edge states as well as the gapped bulk states [3-l3|- In 
order to have topological phase, one ought to have an 
gap in the energy spectrum separating the ground state 
from the excited states. In such a case, one can define 
a smooth deformation in the Hamiltonian which does 
not close the bulk gap. This is due to the fact that one 
gapped state can not be deformed into another gapped 
state in a different topological class unless a quantum 
phase transition occurs when the system become gap¬ 
less. This field has attracted a lot of interest due to 
its wide range of applicability in various areas of con¬ 
densed matter systems such as, quantum Hall effect Q, 
superconductors Impl , Z 2 topological insulators (spin 
Hall insulators) [lll - ll^ etc. The topological phases are 
crucially dependent on some particular symmetries of the 
system such as time reversal(TR), space inversion(SI), 
particle-hole(PH) and chiral etc. The gapless edge states 
are topologically stable against those perturbations that 
do not break the symmetries of the system. The topolog¬ 
ically protected gapless edge states play important role 
in determining the transport properties of the system. 
The total number of topologically protected edge modes 
in a given system is associated with the topological 
numbers such as the Thouless-Kohmoto-Nightingale-den 
Nijs(TKNN) number (the first Chern number) for the 
systems without time reversal symmetry!^ Q and the 
Z 2 invariant in case of time reversal invariant systems [ll|. 

The two-dimensional(2D) topological insulators were 
theoretically predicted by Bernevig et. al. [3 El and ex¬ 
perimentally observed in HgTe/CdTe quantum wells El- 
Such an insulator was already proposed by Kane and 
Mele in 2005 El- These quantum states of matter belong 
to a class which is invariant under TR symmetry and 
the spin-orbit (SO) coupling is essential to achieve this. 
Soon after, it was generalized to superconductors and 
superfluids EMU- In 2D, the classification of topologi¬ 
cal SC is similar to that of topological insulators. For 
example, the TR breaking SC are classified by an integer 


ACfl, El similar to that of quantum Hall insulators d 
whereas TR invariant SC are represented by a Z 2 invari¬ 
ant in 2D and IDEB^II- The TR breaking topological 
SC have attracted a lot of attention recently due to their 
relation to non-Abelian statistics [l,ll and their potential 
application to topological quantum computation!^. 
The nature of the low-energy gapless edge states in such 
systems are non-trivial. They imply fractionlization of 
quasi-particles El as well. For example, in a vortex 
core of a spinless pip SC, the zero mode is described 
by a M^orana fermion which is half of a conventional 
fermion^, El- Such a vortex with a Majoranan fermion 
obeys non-Abelian statistics which is crucial for the 
construction of fault-tolerant quantum computers. The 
existence of zero energy Majorana mode 2^ in a vor¬ 
tex core characterizes the topological order in the system. 


The search for the possible realization of topological 
phases in condensed matter systems is an intriguing and 
challenging issue. This involves novel concepts as well as 
potential applications. In this communication, we con¬ 
sider a coexistence phase of singlet SC and SDW which 
induces a triplet SC component. This is precisely due 
to the non-commutativity of the former two orders. The 
singlet SC is taken to be of dy; 2 _y 2 symmetry whereas 
the SDW order parameter is of s-wave and the triplet 
SC is of p -|- ip type symmetry. Such a Hamiltonian is 
shown to yield a non-trivial coexistence phase which is 
topological in addition to the conventional one. A phase 
diagram characterizing different topological phases is 
construted. The Chern numbers and hence the nature of 
the topological phases are determined. The edge state 
spectrum and the vortex state in such system are also 
discussed. 


I. THEORETICAL FORMULATION 

Motivated from the recent spectroscopic experimen¬ 
tal results on the appearence of a nodal gap on the 
deeply underdoped cuprate scEiil , there has been few 
studies [343| to uncover such new and unexpected re- 
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suits. In an earlier work, we have already discussed about 
the topological study of d-wave SC in presence of SDW 
order and compared with the above cuprates datalSSj. 
However, from the study of the group algebrafsl, l35j|. it 
is known that the coexistence of any two non-commuting 
order parameters produces a third order parameter. In 
case of SDW and d-wave superconductivity, there is a 

I 
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Here, is the bare dispersion, U is the onsite Coulomb 
interaction, are the pairing strengths for d-wave 

and p-wave superconductivity and N is the number 
of sites. We model the bare dispersion in the tight- 
binding approximation on a 2D square lattice as = 
—2t{coskx -\- cosky) — it'cos kj; cosky — p. (cka) 

denotes creation (annihilation) operator of the electron 

with spin a = (t,4-) at k = {kx,ky). Here, J2k 
the sum of k over the reduced Brillouin zone (RBZ). 
We express the wave-vector k in units of ■^, with ’a’ 
the lattice parameter of the underlying square lattice. 
Q = (7r,7r) is the SDW nesting vector in 2D. We as¬ 
sume here a commensurate SDW so that k -1- Q = k — Q. 
The staggered spin magnetization is defined as Mq = 
Hk^a^ < ''k+Q >• Since we will be discussing 

about three order parameters below, the crystal symme- 


third, dynamically generated, order parameter [s^. This 
happens to be a triplet SC here. In the present work, 
we consider the coexistence of SDW and d-wave super¬ 
conductivity which can generate a triplet and non-zero 
center of mass superconducting order parameter. We, 
thus start with a Hamiltonian on a 2D square lattice as. 


'-QqCfe'.t + V^{k, k')cl Y-k,i^-k',iCk',^ 

k,k' 

+ )cl Y-k-Q,i''-k'-Q,i'^k',-f- (1) 

k,k' 


try of them should be such that the commutator of any 
two of them should give the third one. For this rea¬ 
son, if is assumed to be of singlet d-wave symmetry, 
the SDW state guarantees that should be of triplet 
type. So we get the singlet interaction = V^SkSk' 
and = V^pkPk', where Sk = ^{coskx — cosky) and 

Pk = sinkx + i sin ky. We assume that are attrac¬ 
tive. The SC order parameters are defined as, for sin¬ 
glet state, Al, = Ajsfc/ = Vo^Sk''E/,Sk < >= 

other hand, 

the triplet order is decoupled in the main band and in 
the magnetic band as A^, = A^pk' = V^pk' J2k Pk < 
ci+Q^Y-kA >= ^lisin kx+i sin ky) = Af ^-hiA^ ^ while 

^k* = ^oPk' Y^kPk < >■ ^® reason 

why the time-reversal symmetry remains invariant in the 
triplet SC state. After substituting these mean-field or¬ 
ders, the total Hamiltonian reads as. 


n = 


k,(7 


ikc\^^Ck,a + Mo ^-kA^k^) + A/^cl Y-k-QA ^k C-fc-Q.4,Cfc,t 


k,k' 




^k+Q,a^k+Q,a) + ^k (''k,tT'^k,a ^k+Q^a-'^k+Q,a)] + Mq ^a'^k,(y 

k,(T 


+ ^[^/c(Cfc,-|-cLfc_4, + C-feqCfe^'l- c\j^QY-k-QA ^-k-QA^k+QA)] 
k 


+ ^\+QA^^-kA^ ^k {c—k-QA^kA C-fc^Cfe+Qq)],(2) 

k 


where = —it' cos kx cos ky — p and = 

—2t{coskx + cosky). In the above Hamiltonian, 
the nesting property in the band dispersion i.e. 


i+Q = ^k^ 4+q = and the order param- 

COS kxCi—COS ky a \ /\ 2 


eters A^g 


= -Ai( 


= —Al and 


^k+Q ~ ~^k have been employed. In the momen¬ 
tum space, the Hamiltonian can be expressed as, 
ijjl'H{k)ipk where the four-component spinor 
ipk is, V’i = {.cl-\A-k-Qi,cYAk+Q^). Thus, the 
Hamiltonian matrix 1-L{k) in this basis is written as, 
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In what follows, we study the energy spectrum of the 
above Hamiltonian. The Hamiltonian (eqn.(3)) is diago¬ 
nalized and the quasiparticle spectrum is obtained as, 


E±,±ik) = ± 



(Ar^+|A2 p +Mi±2^i 


-2 
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FIG. 1. Energy spectra E±^+{k), corresponding to coexis¬ 
tence of SC order parameters d^ 2 _y 2 -I- (px + ipy) and that 
of the SDW order parameter showing fully gapped spec¬ 
trum. For illustration, here, we have chosen Mo = 0.3t eV, 
t’ = -0.3t, Aj = Ag = t, /i = 0.25t {t = 0.15 eV.) 

It is obvious that the energy spectrum is fully gapped 
and the gap closes only when the following condition is 
satisfied. 


explore the Chern number associated with such phases 
in the next section. 

Ml 



FIG. 2. Schematic phase diagram of the coexistence phase 
of SC d^ 2 _y 2 + {px + ipy) order and SDW order in the case 
where k = (0,0) and (7r,7r). The a:-axis labels the parameter 
H and the y-axis labels the parameter Mq. Integer Af labels 
the Chern number of the coexistence of SC and SDW. 


^+^ + ^-^ + iAir+\Al\^ +M^ 

/- (4) 

=I I' Ai -mtr+ 4^'4“'- 

It is found from the straight forward calculation [s^ that 
this condition is equivalent to 

4+2 + (A1)2+m2=4-2+|A2 p, 

{AlC + Mo\Al 1)2 + (Al)2(| a2 I -Al,r (5) 

+2 Al (Al /. Al —Mo^j!")(| A^ I —Al^/.) = 0. 

The second equation in Eqn.(5) is met only when k = 
(0, 0) and (tt, tt). These two points give the same follow¬ 
ing condition, 

16t^+ M^ = (4t'+pf. (6) 

When the above condition (Eqn.(6)) is satisfied, the en¬ 
ergy gap closes. Thus, we find the topologically trivial 
and the non-trivial regions as shown in Fig.2. We will 


II. PHASE DIAGRAM AND THE CHERN 
NUMBER 

To examine the topological phase transition(TPT), it 
is convenient to use the dual Hamiltonian instead of the 
original one. This can be done through a constant uni¬ 
tary transformation matrix D as, 

H°{k) = DHik)D\ D = E( 'j , ( 7 ) 

resulting 

/ 4 + + M 0 Al-Al 0 C 

y^(k)= A 2 *-A 1 -C-Mo - 4 - 0 

^ ’ 0 - 4 - - 4 + + M 0 -Al*-Al 

V 4” 0 -Al-Al 4 +-Mo 

( 8 ) 

In the limit t —>■ 0, the dual Hamiltonian 'H^(k) in the 
leading order around k = (0,0) and (tt, tt) gives rise to 
the following two 2x2 block-Hamiltonian as, 
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/-At'-fj, + Mo Al{ka; + iky) 
Al{kj:-iky) 4<' + ^-Mo 
0 0 

Vo 0 


We notice here that the above dual Hamiltonian have 
a close similarity to the Hamiltonian of the spinless chiral 
p + ip supercondutor discussed inQ. 

In order to study the phase diagram of this Hamilto¬ 
nian one needs to determine the phase boundaries corre¬ 
sponding to gapless regions since the topological invari¬ 
ants can not change without closing the bulk gap. For the 
present model, the critical lines are determined by solving 
equation(8), i.e., Mq = ±(4t' -|- p) for the upper(lower) 
blocks in the case with k = (0,0) and (7r,7r). The phase 
becomes topological in the region when Mq < ±(4t' -|- p) 
whereas it is trivial for Mq > ±(4t' -|- p). 

Similarly, one can draw a phase diagram for any k 
point in the RBZ. e.g. for k = (0,7r) and (7r,0), we get 
the same condition for the topological phase transition 
to occur from second condition of Eqn.(5), i.e. (4t' — 
p) 7 ^ 0,AJ = 0. This leads to the following Eqn. of 
the critical lines Mq = {At' — pY ■ The topological phase 
exist in the region Mq < ±(4t' — p) whereas it is trivial 
for Mq > ±(4t' — p). 

Based on the finiteness of the Chern number given be¬ 
low, we propose a phase diagram in the in the former 
case (Eig. 2) in the {p,Mq) plane which distinguishes 
the topological and non-topological phases. This is the 
new result of the present manuscript. Thus, it is obvi¬ 
ous that both these phases are separated by a quantum 
phase transition line. The Chern numbers in the topo¬ 
logical phases of Eig. 2 are calculated below. 

It is well known that the topological phases can be 
characterized by Chern numbers. Eor a specific model 
Hamiltonian h{k) = da{k)aa, with ctq, the Pauli ma¬ 
trices and da{k) = [di{k),d 2 {k),d 3 {k)], the Chern num¬ 
ber can be calculated from the expression 
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where the unit vector d{k) = (i{k)/^Y^ d?{k) charac¬ 
terizes a map from the Brillouin zone vector k to unit 
sphere. The Chern number simply counts the number of 
times d{k) wraps around the unit sphere as a function of 
k. In the present model the Chern number for the case 
Mq < ±(4t' -I- p) is calculated as A/” = 1 whereas it van¬ 
ishes when Mq > ±(4t' + p). In presence of SDW order, 
the p + ip SC state has odd parity and Sz=0 symmetry 
which is a fully gapped system. Due to SDW order, it 
has 17(1) spin rotation along say, 2 :-axis and 7ro(C'2) = Z 
which corresponds to class A in accordance with the sym¬ 
metry classification of Altland and Zirnbauerjs^. This 
means that there are infinite number of distinct topologi¬ 
cal SC with Sz conservation and are labeled by an integer 
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FIG. 3. Edge-state spectrum of the coexistence phase of px -I- 
ipy SC in presence of SDW order on a cylindrical geometry. 
Parameters are chosen as, t' = —0.30t eV, Aq = t, p — .25t, 
Mq = .3t for a lattice of W = 100 sites. 


which is the Chern number [slj. This is associated with 
the number of chiral fermion edge modes. 


III. EDGE STATES AND THE VORTEX 
STRUCTURE 

In order to see the evolution of edge states in this 
model(coexistence oip+ip SC and SDW order), we stud¬ 
ied it numerically on a cylindrical geometry with peri¬ 
odic boundary condition in y-direction and open bound¬ 
ary condition in x-direction. We solved the eigen-value 
problem where the Hamiltonian has been diagonalized 
on Nx=l00 sites. The energy dispersion Ek versus ky 
has been obtained and hence the edge states has been 
shown in Fig. 3. As it is already been discussed in the 
previous section, two chiral edge states characterize the 
topological phase in this model. 

It is well known that the vortex of a topological SC 
with odd topological quantum number M carries an odd 
number of Majorana zero modes. The existence of such 
zero modes in the vortex core oi & p + ip SC is shown 
to be due to index theoremjH, In the present case, 
since the SC is coexisting with SDW order, the topo¬ 
logical classification is always trivial for class A in ID. 
Since the existence of the zero modes in the vortex core 
is determined by the symmetry classification in one space 
dimension less, the topological p -I- ip SC with the chiral 
edge modes won’t support such zero energy vortex bound 
stated. 
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IV. CONCLUSION 


In conclusion, we summarize the main findings of the 
present manuscript. We consider a possible coexistence 
of singlet SC and SDW which induces a triplet SC com¬ 
ponent as well. The singlet SC is taken to be of d^ 2 _y 2 
symmetry whereas the SDW order parameter is of s-wave 
and the triplet SC is p + ip type symmetry. Such a 
Hamiltonian is shown to yield a non-trivial coexistence 
phase which is topological in addition to the conventional 


one. A phase diagram characterizing different topological 
phases is construted. The Chern numbers and hence the 
nature of the topological phases are determined. The 
edge state spectrum and the possibility of whether the 
vortex state harbouring the zero modes are discussed. 

V. ACKNOWLEDGEMENTS 

This work is supported by Council of Scientific and 
Industrial Research (CSIR), India. 


[1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ). 

[2] Xiao-Liang Qi and Shou -Cheng Zhang, Rev. Mod. Phys. 
83, 1057 (2011). 

[3] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 

[4] C. Nayak and F. Wilczek, Nucl. Phys. B 479, 529 (1996). 

[5] E. Fradkin, C. Nayak, A. Tsvelik and F. Wilczek, Nucl. 
Phys. B 516, 704 (1998). 

[6] Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993). 

[7] D. J. Thouless, M. Kohmoto, M. P. Nightingale and M. 
den Nijs, Phys. Rev. Lett. 49 405 (1982). 

[8] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

[9] D.A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 

[10] A. Stern, F. von Oppen and E. Mariani, Phys. Rev. B 
70, 205338 (2004). 

[11] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

[12] L. Sheng, D. N. Sheng, C. S. Ting and F. D. M. Haldane, 
Phys. Rev. Lett. 95, 136602 (2005). 

[13] B. A. Bernevig and S. -C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006). 

[14] B. A. Bernevig, T. L. Hughes and S. -C. Zhang, Science 
314, 1757 (2006). 

[15] M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X. -L. Qi, and S. -C. Zhang , 
Science 318 766 (2008). 

[16] R. Roy, arXiv:0803.2868 ]cond-niat,2008]. 

[17] A. P. Schnyder, S. Ryu, A.Furusaki and A. W. W. Lud¬ 
wig, Phys. Rev. B 78, 195125 (2008). 

[18] A. Kitaev, AIP conf. Proc. 1134, 22 (2009). 

[19] X. -L Qi, T. L. Hughes, S. Raghu and S. -C. Zhang, Phys. 
Rev. Lett. 102, 187001 (2009). 

[20] G. E. Volovik, Zh. Eksp. Theor. Fiz. 94, 123 (1988) [Sov. 
Phys. JETP 67, 1804 (1988)]. 

[21] C. Nayak, S. H. Simon, A. Stern, M. Freedman and S. 
Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

[22] D. H. Lee, G. M. Zhang and T. Xiang, Phys. Rev. Lett. 
99, 196805 (2007). 

[23] N. B. Kopnin and M. M. Salomma, Phys. Rev. B 44, 
9667 (1991). 

[24] G. E. Volovik, Pis’ma Zh. Eksp. Theor. Fiz. 70, 601 
(1999) [JETP Lett. 70, 609 (1999)]. 


[25] K. Tanaka, W. S. Lee, D. H. Lu, A. Fujimori, T. Fujji, 
Risdiana, I. Terasaki, D. J. Scalapino, T. P. Devereaux, 
Z. Hussain, Science 314, 1910 (2006). 

[26] I. M. Vishik, M. Hashimoto, R. -H. He, W. -S. Lee , F. 
Schmitt, D. Lu, R. G. Moore, G. Zhang, W. Meevasana, 
T. Sasagawa, Proceedings of the National Academy of 
Sciences 109, 18332 (2012). 

[27] A. Ino, G. Kim, M. Nakamura, T. Yoshida, T. Mizokawa, 
Z. X. Shen, A. Fujimori, T. Kakeshita, H. Eisaki and S. 
Uchida, Phys. Rev. B 62, 4137 (2000). 

[28] E. Razzoli, G. Drachuck, A. Keren, M. Radovic, N. C. 
Plumb, J. Ghang, Y. B. Huang, H. Ding, J. Mesot and 
M. Shi, Phys. Rev. Lett. 110, 047004 (2013). 

[29] Y. Peng, J. Meng, D. Mou, J. He, L. Zhao, Y. Wu, G. 
Liu, X. Dong, S. He, J. Zhang, Nature Communications 
4, 2459 (2013). 

[30] K. M. Shen, T. Yoshida, D. H. Lu, F. Ronning, N. P. 
Armitage, W. S. Lee, X. J. Zhou, A. Damascelli, D. L. 
Feng, N. J. C. Ingle, H. Eisaki, Y. kohsaka, H. Takagi, T. 
Kakeshita, S. Uchida, P. K. Mang, M. Greven, Y. Onose, 
Y. Taguchi, Y. Tokura, S. Komiya, Y. Ando, M. Azuma, 
M. Takano, A. Fujimori, Z. X. Shen, Phys. Rev. B 69, 
054503 (2004). 

[31] Yuan-Ming Lu, Tao Xiang and Dung-Hai Lee, Nature 
Physics 10, 634 (2014). 

[32] Tanmoy Das, arXiv.1312.0544 vl ]cond-mat.supr-con] 02 
Dec. (2013). 

[33] Amit Gupta and Debanand Sa, Solid State Gommun. 
203, 41 (2015). 

[34] S. -G. Zhang, Science 275, 1089 (1997). 

[35] R. S. Markiewicz and M. T. Vaughn, Phys. Rev. B 57, 
R14052 (1998). 

[36] B. Kyung, Phys. Rev. B 62, 9083 (2000). 

[37] M. Sato, Phys. Rev. B 73, 214502 (2006). 

[38] A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 
(1997). 

[39] Sumanta Tewari, S. Das Sarma and Dung-Hai Lee, Phys. 
Rev. Lett. 99, 037001 (2007). 

[40] Masatoshi Sato and Satoshi Fujimoto, Phys. Rev. B 79, 
094504 (2009). 



